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Abstract. In this paper, we review tests of the strong equivalence principle (SEP) 
derived from pulsar-white dwarf data. The extreme difference in binding energy 
between both components and the precise measurement of the orbital motion provided 
by pulsar timing allow the only current precision SEP tests for strongly self-gravitating 
bodies. We start by highlighting why such tests are conceptually important. We then 
review previous work where limits on SEP violation are obtained with an ensemble 
of wide binary systems with small eccentricity orbits. Then we propose a new SEP 
violation test based on the measurement of the variation of the orbital eccentricity 
(e). This new method has the following advantages: a) unlike previous methods it is 
not based on probabilistic considerations, b) it can make a direct detection of SEP 
violation, c) the measurement of e is not contaminated by any known external effects, 
which implies that this SEP test is only restricted by the measurement precision of 
e. In the final part of the review, we conceptually compare the SEP test with the 
test for dipolar radiation damping, a phenomenon closely related to SEP violation, 
and speculate on future prospects by new types of tests in globular clusters and future 
triple systems. 



1. Introduction 

The Strong Equivalence Principle (SEP) extends the Weak Equivalence Principle (WEP) 
to the universality of free fall (UFF) of self-gravitating bodies. General Relativity (GR) 
assumes the validity of WEP and SEP, i.e. in GR the world line of a body is independent 
of its chemical composition and gravitational binding energy. Therefore, a detection of 
SEP violation would falsify GR. On the other hand, alternative theories of gravity 
generally violate SEP. This is also the case for most alternative metric theories, such 
as the Jordan-Fierz-Brans-Dicke (JFBD) theory [321 E31 [7] or the more general scalar- 
tensor theories of gravity [11], which satisfy WEP as a consequence of their postulate 
of universal coupling between matter and gravity ([61], see also review by T. Damour 
in this issue j3- 

| In fact,Willy Scherrer was the first onewho had proposed a scalar-tensor theory of gravitation (see 
[29] for details on the history of scalar-tensor theories). 
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SEP violation could have observable consequences even in the weak, quasi- 
stationary gravitational fields in our Solar System [391 HQ], i n particular a "polarization" 
of the Earth-Moon system by the external Solar field [41] . Detecting such a polarization 
is one of the main motivations of the Lunar-Laser- Ranging (LLR) experiment, which is 
described in other articles in this issue. 

In view of the smallness of self-gravity of Solar System bodies, the LLR experiment 
says nothing about strong-field aspects of gravity, like the gravitational properties of 
extremely compact objects like neutron stars or black holes. In the presence of such 
bodies, alternative theories of gravity that pass Solar System tests can still produce 
various observable phenomena not predicted by GR, among these the radiation of dipolar 
gravitational waves and strong-field SEP violation (Section |2J). The former has not yet 
been detected in the most sensitive binary pulsar experiments [371 EE] , which given the 
current experimental precision already provide extremely stringent tests of alternative 
theories of gravity. In the near future the detailed study of the gravitational wave 
signal emitted by mergers of compact objects will also provide independent tests of the 
radiative properties of gravity (e.g., [62]). 

For strong-field SEP violation, the best current limits come from millisecond pulsar 
(MSP) - white dwarf (WD) systems with wide orbits. If there is a violation of UFF by 
neutron stars, then the gravitational field of the Milky Way would polarize the binary 
orbit [16] . Such tests of the UFF of strongly self-gravitating masses are the subject of 
this work. In comparison with the LLR tests, they have two disadvantages, one of them 
being the much weaker polarizing external field (|g| ~ 2 x 10 _10 ms -2 compared with 
the acceleration of the Solar gravitational field at the Earth, ~ 6 x 10 -3 ms -2 ) and the 
precision of the ranging, which is of the order of 10 m for the best pulsar experiments 
(10 _2 m for LLR). This is almost completely compensated by the gravitational binding 
energy of neutron stars, -Eg rav) which is a large fraction of the total inertial mass-energy: 
£ g rav = E gTav /Mic 2 ~ —0.15; this is more than eight orders of magnitude larger than 
Earth's (e gr av,e ~ — 5 x 10~ 10 ). This results in experiments with comparable limits on 
SEP violation, which are nonetheless complementary since they probe different regimes 
of binding energy. 

In Section [2] we summarize some of the theoretical foundations of SEP violation for 
strongly self-gravitating bodies. In Section [3] we review the D amour- S chafer test, which 
yields the best current limits on SEP violation by strongly self-gravitating bodies, if 
applied to a whole population of small-eccentricity systems. 

In Section H] we suggest new pulsar timing experiments that avoid the probabilistic 
considerations of present tests and have the potential to detect SEP violation; these 
attempt to directly measure its effects, in particular the variation of the orbital 
eccentricity, e. It is shown that this measurement is not contaminated by external effects 
and, because of this, the limits on SEP violation are only restricted by the precision 
of the measurement of e. Our simulations show that such direct tests will very soon 
surpass the best current limits on SEP violation for strongly gravitating bodies. 

Section [5] provides a short discussion on the complementarity of SEP tests and tests 
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for dipolar radiation damping in constraining alternative theories of gravity. Finally, in 
Section [6] we briefly discuss the possibilities of the SEP test proposed Section 0] when 
the test binaries are accelerating in a fields much stronger than that of the Galaxy, 
like that of a globular cluster or an additional outer companion. Such systems would 
likely offer very significant gains in the power of SEP tests, essentially providing an 
experiment with fields and velocities very similar to those present LLR experiment, but 
with "proof masses" with a ~ 10 8 larger binding energy - essentially a combination of 
the best features of both types of experiments. 



2. SEP violation and strongly self-gravitating bodies 

2.1. Beyond the weak- field approximation 

In general, an alternative theory of gravity is expected to violate the strong equivalence 
principle, in the sense that the ratio between the gravitational mass Mq of a self- 
gravitating body to its inertial mass M\ will admit an expansion of the type 

Mq 2 

■j^ = 1+A= 1 + 7?i Sgrav + 7? 2 £ grav + . . . (1) 

At the first post-Newtonian level r\\ parametrizes the weak-field violation of SEP 
(Nordtvedt effect). In the Parameterized Post-Newtonian (PPN) framework this 
Nordtvedt parameter is given as a combination of different PPN parameters (see [61] 
for details on the PPN formalism): 

10 2 2 1 

T)! = T) N = 4p - 7 - 3 - -j£ - «i + ~a 2 ~ -Cl - gC2 • (2) 

The parameter t]n is well constrained by lunar laser ranging (LLR) experiments in 
the Solar System (see contributions on LLR in this review). In view of the fact that 
£ g rav,e ~ -5 x 1CT 10 , e grav ,Moon ~ -2 x 10~ u and £g r av,0 ~ -1CT 6 , it is clear that 
higher order deviations from SEP cannot be tested in LLR experiments, or any other 
experiment in the Solar System, in the foreseeable future. To test a violation of SEP that 
might occur beyond the weak-field regime, one needs strongly self-gravitating bodies. 
Presently, in nature the best objects for such tests are neutron stars observed as active 
radio pulsars (e grav ~ — 0.15) j§| 

Since beyond the weak field approximation there is no general PPN formalism 
available, discussions of gravity tests in this regime are done in various theory specific 
frameworks. An example for a very detailed investigation of higher order /strong-fied 
deviations from GR, within the family of (well-defined) scalar-tensor theories of gravity, 
are the frameworks developed by Damour and Esposito-Farese. Specifically, 

A) The 4-parameter framework T (j, (3;e,() of [IB], which defines the second post- 
Newtonian extension of the original (Eddington) PPN framework To (7, 0). In this 



§ For gravity theories that predict a no-hair theorem for black holes, neutron stars are in fact the most 
compact objects suitable for strong-field SEP violation tests (see e.g. [131 S21 D] ) • 
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framework 

A* ~ -~(4/3 - 7 - S)c A + (| + C) b A , (3) 

where the compactness ca ~ — 2e gra v, J 4 and ~ c A . 

B) The 2-parameter class of bi-scalar-tensor theories T2((3',/3") introduced in pTj . 
Within this framework, one has no Nordtvedt effect in the Solar System, since 
i] N = 0. On the second post-Newtonian level one has e = (3', ( = 0. And (3" 
parametrizes contributions beyond the second post-Newtonian level. 

C) The 2-parameter class of mono-scalar-tensor theories T\(ao, (3q) of Q3]> which 
for certain values of (3o exhibits significant strong field deviations from GR, and a 
corresponding violation of SEP for neutron stars. 

For details we refer the reader to [TO] and references therein. In the following we will have 
a closer look at case C, as this class of scalar-tensor theories of gravity illustrates quite 
impressively a violation of SEP that can only be measured with strongly self-gravitating 
bodies. We would like to note, however, that this particular case of scalar-tensor 
theories of gravity resembles just one example of how non-linearities in the gravitational 
interaction could drive gravity away from GR in the strong fields of compact masses 
and lead to a strong-field violation of SEP. 

2.2. Strong field effects and the violation of SEP. 

Damour and Esposito-Farese [12] found that scalar-tensor theories, which pass the 
weak-field tests in the Solar System, could still exhibit large, strong-field-induced 
deviations in systems involving neutron stars ("spontaneous scalarization"). This has 
been studied extensively in a two- parameter space of theories, Ti(a ,/?o), defined by 
the coupling function which is a quadratic polynomial in the scalar field ip: a(ip) = 
ao(ip — ipo) + (3o((p — V 9 o) 2 /2 p2J ttH US]- The parameter a defines the linear matter- 
scalar coupling constant and (3 the quadratic coupling of matter to two scalar particles, 
while higher-order vertices are neglected. In this sense, this is a natural extension of 
JFBD gravity. 

In presence of such non-perturbative strong-field deviations away from GR, we can 
have a situation where the effective coupling strength of the neutron star, a a, is of order 
unity, even if the scalar- matter coupling, ato, is unobservably small in the Solar Syste n§ 
Such an effect leads to a violation of SEP that requires test systems which contain a 
neutron star. 

The structure dependence of the effective gravitational constant Gab, has the 
consequence that the pulsar does not fall in the same way as its companion in the 

|| The quantity a a = dlnMA/dtpo measures the effective strength of the coupling between a self- 
gravitating body A, with total mass Ma, and the scalar field ip. It is equivalent to the negative ratio 
of total scalar charge to total mass. For a weakly self-gravitating body a a — «o- 
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gravitational field of a third body, which in our case is the Galaxy. One finds for a 
pulsar with a weakly self-gravitating companion, since oiq -C 1, that [11] 

A p - A c ~ a (a p - a c ) ~ a (a p - a ) . (4) 

While |a | < 0.003 by the Cassini experiment [3], a p can be of order unity for neutron 
stars, as outlined above. Although the effect is greatly suppressed by a small factor a , 
it could in the presence of non-perturbative strong-field deviations be many orders of 
magnitude stronger than in the Solar System (factor a p /a ), and therefore still become 
visible in binary pulsar timing experiments. 



2.3. SEP violation and orbital motion of binary pulsars 

In the previous section we have seen that one could expect violations of SEP in gravity 
regimes where the Solar System tests are insensitive to, and therefore one has to utilize 
test systems that contain strongly self-gravitating bodies. Currently, binary pulsars 
are the best probes for testing such kind of strong field gravity effects [10J. Before we 
discuss in details the SEP tests that can be conducted in timing binary pulsars, we need 
to outline the orbital dynamics of a binary system that falls freely in the gravitational 
field of our Galaxy in the presence of a SEP violation. 

In case of a violation of SEP the equations for the relative motion R(t) are given 

by 

R = -GM^- + A PN + A A (5) 

where M is the total mass of the binary and Q = G pc is the effective gravitational 
constant between the two bodies [16J. The post-Newtonian contributions are denoted 



by A PN and the additional acceleration caused by SEP violation, A a, is given to leading 
order bjjjjl 

A A = (A p - A c ) g , (6) 

where A, is defined by the ratio between "gravitational" and "inertial" mass as given in 
eq. ([1]), and g is the external acceleration caused by the gravitational field of the Galaxy. 
The secular changes to the binary system with orbital frequency n\> and eccentricity e 
caused by a violation of SEP are [16] 

(dn h /dt) = , (de/dt) = f x 1 + cj PN k x e , (dl/dt) = f x e , (7) 

where 

e^ea, 1 = VT^k , f = (A p - A c ) g , (8) 

and Vo = {QMn^) 1 ^ is a measure for the relative orbital velocity between the pulsar 
and its companion. The unit vector a points towards periastron and the unit vector k is 

As the gravitational field of the Galaxy at the location of the binary pulsars is weak, and since most of 
the mass of the Galaxy is made of non strongly-self-gravitating bodies, any higher order contributions 
to A a are negligible. 
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Figure 1. Time evolution of the observed orbital eccentricity vector e = epN + eA 
in a small eccentricity binary. The vector gj_ represents the projection of the external 
acceleration in the orbital plane. 



parallel to the orbital angular momentum. The post-Newtonian precession of periastron 
is given by 

^ (Vo/c) 2 

w PN = 3J 7 '- n h , (9) 

1 — e 1 

where J 7 is a theory-dependent factor [T = 1 in GR). As a result of eqs. ([7J), a violation 
of the UFF would lead to a change in the observed orbital eccentricity 



e = (de/dt) ■ a = Vl - e 2 (b • f) , (b = k x a) (10) 
and a change in the angle i between the line-of-sight direction to the pulsar, Ko, and k 
d cos i d ( K • 1 \ e 



dt dt {y/TT^J y/T 



,2 



(K -b)(k-f). (11) 



In binary pulsar timing experiments this change in i becomes apparent as a change in 
the timing parameter x, which is the projected semi- major axis of the pulsar orbit given 
by x = a p sin i/c. 

As shown in [16], for very small eccentricities the equations of motion ([7]) essentially 
decouple. As a consequence, the orbital plane remains fixed and the evolution of the 
eccentricity can be written as 

e(t) = e PN (t) + e A , e A = f±/wp N . (12) 

The vector epN(t) has a fixed length and is turning in the orbital plane with angular 
velocity ojp^. The polarization of the orbit due to the violation of SEP is represented 
by the constant eccentricity vector eA, which points into the direction of the projection 
of f into the orbital plane (denoted by fj_). Figure Q] illustrates the time evolution of the 
orbital eccentricity. Thus, in the presence of a SEP violation the observed eccentricity 
oscillates between a minimum (|epN — €,a\) and maximum (|epN + ca|) eccentricity with 
a period of 27t/o;pn. 
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3. The Damour-Schafer test 

In 1991, when Damour and Schafer first wrote their paper describing the orbital 
dynamics of a binary pulsar under the influence of SEP violation, only four binary 
pulsars were known in the Galactic disk. Two of these (PSR B1913+16 and 
PSR B1957+20) were clearly inadequate for that test - - not only because of the 
compactness of their orbits, but also because one of them is a double neutron star 
system that lacks the required amount of asymmetry in the binding energy, necessary 
for a stringent test of SEP violation [the term (a p — a c ) in eq. (jlj)]. The remaining 
systems were PSR B1855+09 [15] and PSR B1953+29 ©. The latter was already 
deemed to be the best probe for the detection of SEP violation effects because of its 
orbital period of 117 days, by far the largest then known. Nevertheless, even for this 
system, the orbital effects predicted for allowed levels of SEP violation — in particular, 
the changing eccentricity e — were too small for detection given the timing precision 
and baseline of that system in 1991. Another way of saying this is that at that time no 
interesting limits on SEP violation could be derived from existing constraints of e. It is 
for this reason that Damour and Schafer proposed a novel, indirect statistical SEP test, 
by assuming that SEP violation is responsible for a maximal part of the observed small 
orbital eccentricities. 

3.1. Small- eccentricity binary pulsars and SEP 

The basic idea presented by Damour and Schafer has been already described above, 
namely that the violation of the SEP introduces a polarization of the orbit of binary 
pulsars that is best represented by a vector addition where the observed eccentricity 
vector e(t) lies on a circle (see Fig. [1]). In the extreme case, however, the vector 
addition leads to a near or even complete cancellation of the eccentricity vector and 
a SEP violation is therefore not detectable. Unfortunately, the intrinsic vector epN, 
and therefore also its orientation relative to the SEP component vector eA, is unknown. 
However, for a sufficient age of the system one can assume that the relativistic precession 
of the orbit will have caused the eccentricity vector to have made many turns since the 
system's birth, thereby effectively randomizing the relative orientation 9 = u-p-^t. In fact, 
the angular velocity of the periastron advance, cjpN, should be appreciably larger than 
the angular velocity of the rotation of the Galaxy with which g rotates in the reference 
frame of the binary system. As a result, the projection of the Galactic acceleration 
vector onto the orbit can be considered constant, and statistical arguments based on 
the exclusion of small 6, i.e. possible near cancellation, can be applied. In the ideal 
case, the masses of the system components, the inclination of the orbit and the distance 
to the pulsar are known. The angle describing the orientation of the orbit around the 
line-of-sight, Q, is generally not observable and has to be treated as an independent 
random variable uniformly distributed between and 2tt. 

As the treatment above has been derived for small-eccentricity, long-orbital period 
pulsars, the figure-of-merit for suitable systems considered by Damour and Schafer was 
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large values of P^/e — the large values of P& decrease the relative orbital velocities 
Vo, which increases the amplitude of f (eq. [8]) and the predicted change of eccentricity 
(eq. [7]); the small value of e implies that, in a statistical sense the orbit has been little 
changed by such SEP violation effects. 

They applied their method to the two systems with (by far) the best P^je at the 
time, PSR B1855+09 and PSR B1953+29. While the measurement of a Shapiro delay 
[4"9[ 127] for PSR B1855+09 [IH [31] were available, providing constraints on sin i and 
the component masses, evolutionary arguments were used to constrain the parameters 
for PSR B1953+29. Overall, Damour and Schafer derived 90% confidence level limits 
of |A P | < 5.6 x 1CT 2 and |A P | < 1.1 x 10~ 2 , respectively. 

Six years later, [59] presented an updated analysis, applying this Damour- Schafer 
test to eight systems with large values for P^/e, this excluded SEP violation at a level 
of 5 x 10~ 3 with more than 95% confidence. 

3.2. A population of small- eccentricity binary pulsars 

The analysis presented in [59J has a caveat, by selecting only those systems with the 
best figure-of-merit P^/e. By this, one introduces a selection bias, as it is possible that 
on the one hand a significant eccentricity (which would reduce the figure-of-merit and 
a possible weight in the analysis) is actually the result of an SEP violation, and on the 
other hand the small eccentricities selected are those where by chance 9 is small. (This 
was pointed out to one of us (NW) by Kenneth Nordtvedt.) 

In order to take this into account, [60] presented an updated analysis that included 
all relevant small-eccentricity binary pulsars at that time. Based on extensive Monte- 
Carlo simulations, a large set of simulated (cumulative) distributions of eccentricities 
was compared to the distribution of eccentricities observed in the population of small- 
eccentricity binary pulsars. In this analysis, unknown angles like 9 and Q were 
distributed uniformly between and 2ir. Based on a Kuiper's test like criterion, the 
number of simulated distributions which are in agreement with the observed distribution 
was determined. As a result, [60] found a 95%-confidence limit of |A P | < 9 x 10 -3 . This 
limit is weaker than the one in [59], but certainly more representative by being able to 
consistently include also systems with seemingly worse figure-of-merits. 

In the same spirit, [52] and later [30] presented an updated analysis, using even 
larger samples of pulsar- WD systems, i.e. the population of all known systems that 
are thought to have evolved with similar extended accretion periods. For deriving 
the median-likelihood value of |A P | for each pulsar, they use a Bayesian analysis 
marginalising over the parameters with similar assumptions to those in [TBI [59] . In 
systems like PSRs 0437—4715 and J1713+0747 where the orientation of the orbit can 
be measured [551 EI] , posterior probability density functions were derived appropriately. 
From all available information they obtained |A p | < 5.6 x 10~ 3 (95% C.L.). The 
best present limit was obtained by [30] using 27 binary systems (including additional 
astrometric and mass information), |A p | < 4.6 x 10~ 3 (95% C.L.). 
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4. Detecting SEP violation directly 

Many new advances in pulsar astronomy have occurred in the two decades elapsed since 
the publication of [16]. Three of them are especially important in this context: 

A) Continuing pulsar surveys have discovered many more binary systems, including 
several with significantly wider orbits than PSR B1953+29, i.e., systems more 
suitable for the detection of SEP violation (Section EJ), 

B) Advances in the sensitivity and bandwidth of radio receivers and improvements in 
the time resolution of instrumentation mean that we can time the binary pulsars 
in these systems much more precisely than possible before, 

C) The sheer amount of elapsed time provides much larger timing baselines. 

All of these developments increase the sensitivity to the direct effects of SEP violation, in 
particular e (Section 14. II) . to the point that it should now be possible to derive even more 
stringent limits on SEP violation from them (Section 14. 2p . This has several advantages 
over the statistical method described in Section 

(i) We can actually detect a potential SEP violation by measuring a non-zero e, 
particularly if the same phenomenon is observed for several binary MSPs. With 
the statistical method reviewed in Section [3] we can only estimate upper limits for 
the effect. 

(ii) For wide pulsar-WD systems, only this hypothetical SEP violation can cause a 
measurable e: orbital circularization due to emission of gravitational waves or 
aberration effects cause a change in eccentricity that is many orders of magnitude 
below the current experimental precision on e; furthermore there are no effects due 
to mass loss in the system or its motion (Section 14 .3p . In this sense, this is a clean 
test of the validity of GR. 

(iii) The clean nature of the e test implies that the limit on SEP violation will improve 
at the same rate as the precision of the measurement of e. Thus, if the latter 
improves continuously with time (Section 14. ip . the same will happen with limits on 
SEP violation. 

(iv) Since the test can be done with a single binary system, we do not need to assume 
that A p is the same for all systems, as assumed in the statistical test of Section l3~2l 
Indeed, as pointed out in [TP] , alternative theories of gravity predict that A p is a 
function of the pulsar mass. Therefore, a rigorous analysis requires an accurate 
knowledge of the masses of the pulsars, which for many of the pulsar-WD systems 
are not available. In view of strong-field effects like "spontaneous scalarization" , 
the combination of binary pulsars in a generic SEP test could even be rendered 
meaningless. 

(v) We do not need to restrict our sample to systems with small eccentricities 
(Section S3]). 

We now discuss which of the known binary pulsars are most suitable for this particular 
test. 
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4-1. Figure of merit for detection of SEP violation 

Our simulations indicate that for the parameters affected by SEP violation the 
uncertainty provided by timing is given by: 

5e~8.0x — -I — , (13) 

5x ~ 5.3 x 6t , (14) 

where N is the average number of TOAs per unit time (which we assumed to be 
uniform in our simulations), St is the rms of the TOA residuals and T is the observing 
baseline. This also implies, in general, that 5e ~ 1.5 5x/x. These expressions assume 
the absence of red timing noise. If those effects are present at measurable levels, then 
the improvement with time is slower [36J. 

Using the expressions in Section [2] and the expression for the total mass M derived 
from the mass function (eq. 3.15 in [18] ) we can re- write e as: 

e = §(A p -A c )(b.g)GA-) 1/2 , < 15 > 

where p = xc(l — e 2 )/ sin 2 is the length of the semi-latus rectum of the pulsar's orbit 
and b is the unit vector along its length. With the two last equations we can calculate, 
in the eventuality of SEP violation, the significance of its measurement: 

f * 1.7 x 10- 5 s- 3 / 2 (A p , 3 - A C)3 ) (b ■ g 10 ) % * dX {1 , e . 
be \y ot 1 m c smt 

where m c = M c /M & , gio = g/10 _10 ms -2 , Aj^ = Aj/10 -3 , T y is the timing baseline in 
years, Nd is the number of TOAs per day and St^ is the TOA r.m.s. in fis. Furthermore, 
limits on viable alternative theories of gravity from pulsar- WD systems [5], [28] imply that 
we can assume G/Q ~ 1, where G is Newton's gravitational constant. If there is SEP 
violation, then given enough time T it will eventually be detected to high significance. 



4-2. Current precision of direct SEP violation test 

To evaluate eq. f[T6|) . we must determine (m c sini) with some degree of precision; this 
requires a significant measurement of the Shapiro delay. Apart from this, u and Q 
(the longitude of periastron and position angle of the line of nodes) are necessary for 
determining the absolute orientation of the orbital plane in space. Furthermore, the 
parallax of the system n x is necessary for determining its distance d and its location in 
the Galaxy, which are necessary to estimate the Galactic acceleration g at that location 
and its projection (b-g) along the direction of the semi-latus rectum. With the exception 
of u) all of these parameters require high timing precision. There are only two systems 
for which all these parameters have been measured, PSR J0437— 4715 and J1713+0747. 
The latter has a significantly larger x, which according to eq. (f!6l) increases its sensitivity 
to SEP violation. Thus in what follows we discuss this system in more detail. 
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Pulsar timing can in principle be used to detect very low-frequency gravitational 
waves [46] [21]. Presently, several large-scale projects are attempting to achieve this with 
precise, sustained timing of several MSPs (a pulsar timing array, or PTA, [56], |22j IT9]). 
PSR J1713+0747 is part of all these efforts because it is detectable by all major radio 
telescopes and is currently one of the three pulsars with the smallest St known, 0.1 /xs 
[19J. It has also been precisely timed for two decades, making its T unusually large. 
All of this makes this system extremely sensitive to e caused by SEP violation. Using a 
PTA pulsar for this experiment also means that the SEP violation test proposed here 
demands no extra time allocation — the scientific results can be obtained for free as an 
added benefit of ongoing efforts. 

We now use TOA simulations to estimate what limits on e can be achieved for 
this pulsar. We do this partly because no values and uncertainties for e have been 
published to date for this pulsar (an illustration of how unexpected a measurable e is, 
see Section 1473]) . but also to estimate future limits on its precision. In our simulations, 
we use TOA datasets with uncertainties, number, start and end times given by Table 
1 of [51J and Table 2 of [19]; the latter dataset appears to have started at the end of 
2004 and continues to the present and it completely dominates our simulated dataset. 
For the first 6 years 2368 TOAs had been taken, at an average rate of 1 TOA per day, 
we assume in our simulations that 30 TOAs were obtained in a single session every 30 
days; furthermore we assume that this uniform rate continues at present. 

With these assumptions, we obtain 5e = 1.6 x 10~ 18 s _1 as a realistic uncertainty 
at the end of 2012. If the measured value for e is consistent with zero, then given the 
location and orbital orientation of this binary, this would result in A p < 2 x 10~ 3 (95% 
C. L.), which is already twice as constraining as the current limits from the Damour- 
Schafer test. If we improve the timing precision by a factor of 2 and keep the same 
timing strategy, then we would have to wait until 2030 for the precision of the test to 
increase by one order of magnitude. Significantly faster progress might be achievable 
for MSPs with larger x, or with improvements in timing precision that will be provided 
by the Square Kilometre Array (SKA, [47J). 

4-3. Cleanest binary pulsar experiment 

Using the results in [18], in particular their eq. (2.4c), we see that the orbital eccentricity 
e, being adimensional, is not affected by the Doppler shift D that must necessarily occur 
during the coordinate transformation from the reference frame of the binary to the Solar 
System Barycenter (SSB). Another example of such an adimensional quantity is the 
Shapiro delay parameter s (or q, [27J). Therefore^ follows that e is not affected even 
when there is a change in that Doppler factor, Zo At any moment the observed e is 
the "intrinsic" e, plus a small term due to geodetic precession. 

It follows from this that if the observed e is too small for detection, the same limit 
will also apply to the intrinsic e, irrespective of D. This is extremely important: as an 

+ D and D arise from the first and second derivatives of the distance to the system, d [50] . 
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example, the intrinsic variation of the orbital period (-P b int ) or the variation of projected 
semi-major axis of the orbit (x int ), are always "polluted" by D [17]. and we thus observe 





D< VJ ={-J -D< (17) 

thus uncertainties in the estimation of D impose a fundamental limit on the precision of 
the measurement of P b mt and therefore ultimately limit the precision of radiative tests 
of gravity [58] . 

Another effect that changes P b Int is mass loss from the binary, which has a lower 
limit given by the loss of rotational energy from the pulsar [T7J ; furthermore outgassing 
from the companion can in some binaries be so large that no test of GR can be made. 
However, such instances of steady mass loss do not affect the orbital eccentricity [31]. 

For very wide MSP- WD systems, tidal effects — in particular orbital circularization 
- are also expected to be extremely small, since both objects are extremely small 
compared to the size of the orbit. They behave effectively as two point masses. 

Finally, if GR is the correct theory of gravity, then e is given by [18J: 

e GR = e GW + e A , (18) 

where the first term is caused by gravitational wave emission and the second is caused 
by a change in the aberration parameter that is to be expected from geodetic precession. 
The leading term of e GW is given by [42] : 



GW _ _ 304 8/3 5/3 g e(l + 121eV304) 

6 " 15 b [g c) (g+l)V3 (l- e 2)5/2 > ^> 

where q = m p /m c is the mass ratio and T & = GM Q /c 3 = 4.925490947 fis is one Solar 
mass in time units. We derive the second term from the equations in [IB] : 

e A T ° m < Jf^ TJT^v^' A ' »?) ( 20 ) 

J(i, A, rf) = — yt (si n i cos ^ sin2r/ + cos i sin A cosr^) , (21) 
sin A 

where 7] is the longitude of the projection of the pulsar spin axis in the plane of the sky 
measured from ascending node and A is the angle between the pulsar spin axis and the 
line of sight from the pulsar to the Earth (this cannot be zero, otherwise there would 
be no pulsations). 

For PSR J1713+0747, these terms are given, respectively, by e GW = —8.2 x 
10~ 29 s -1 and e A = 8.2 x 10~ 28 s _1 J(i, A, rj). Since J(i,\,r)) must be of the order 
of unityEI this means that e GR is at least 9 orders of magnitude smaller than the current 
be. Given the smallness of all polluting terms, we reach the conclusion that the limits 
on A p will improve as much as the experimental precision of e. 

* This is a fully recycled pulsar, therefore the vast majority of its rotational angular momentum came 
from orbiting material. This implies that the spin and orbital angular momenta must be very closely 
aligned, therefore r\ ~ n/2 and J(\,r)) ~ (i.e., much smaller than unity). This means that for this 
kind of system the only likely contribution to e comes from e , which is ten orders of magnitude 
smaller than the present Se. 
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Figure 2. Limits on A p (95% C. L.) as a function of the longitude of the ascending 
node, ri, for the PSR J1903+0327 binary system, using the values for x and e in [2"5] 
and assuming an orbital inclination of 77.47° (the other possible inclination is 102.53°, 
for which we would have similar constraints, but at different values of f2). A distance 
of 6.4 kpc was assumed. The limits in blue are derived from the measurement of x 
and the limits in black from the measurement of e; the red line highlights the overall 
upper limit. 

4-4- Testing SEP with the 1903+0327 binary system. 

An important advantage of these direct tests of SEP violation is that we do not 
need to restrict our study to binaries with small eccentricities. The binary system 
PSR J1903+0327 [H [25] has good timing precision (St ~ 1 /is) and a wide (x = 
105.593 lt-s) orbit for which e = (14 ± 6) x lCr^s" 1 and i obs = +21(3) x lfr 15 lt- 
ss -1 have been published [25J. A precise distance is not known yet, nor Q. However, 
we can use the parameters we know precisely (position in the sky, sinz, uS) and a model 
of the Galactic potential to estimate limits on A p for any assumed d and Q. These are 
displayed graphically in Fig. [2] for d = 6.4 kpc. 

Note that we can estimate A limits for every value of Q because for such eccentric 
systems SEP violation also causes a torque in the orbital plane (eq. ITT]) . This would 
add a component xsep to the kinematic contribution expected for the system at each 
particular Q, x kin [2j [35] and the Doppler contribution for each particular d (eq. ITT]) . such 
that xsep — ^obs — ikin + xD/D. These limits are not as constraining as those derived 
from PSR J1713+0747, but they are already interesting because they apply to a massive 
neutron star (M p = 1.667 ± 0.021 M , 99.7% C.L.), for which no precise constraints of 
A p have been obtained until now (we elaborate on this in Section [5]). Furthermore, this 
limit will increase fast because a) the e and x limits were obtained with only two years 
of data (from 2008 to 2010), so significant improvements will be possible in a relatively 
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short timescale and b) new broadband coherent dedispersion systems will yield a major 
improvement in the timing precision of this system. 

Future infrared interferometric experiments like GRAVITY |3J will be able 
to measure precisely the astrometric motion of the main sequence companion to 
PSR J1903+0327, allowing therefore a precise, independent measurement of ir x , i and 
Q, thus allowing for unambiguous SEP tests in this system. Moreover, if indeed SEP is 
violated at a measurable level in PSR J1903+0327, the eccentric nature of this system 
allows a unique cross-check of this since both xsep and e should be influenced in a 
characteristic way: According to eqs. ( TTUj) and (fTT|) . the quantity isEp/(xe) only depends 
on the orientation and the eccentricity of the pulsar orbit. 

5. The complementarity of SEP violation and dipolar radiation tests 

As discussed in the Introduction, a theory of gravity which predicts a violation of SEP 
is also expected to predict the emission of dipolar gravitational radiation in asymmetric 
binary systems, like pulsar- WD binaries [61]. By now there are several binary pulsars 
that provide tight constraints on the existence of dipolar gravitational radiation within 
scalar-tensor gravity, as well as within more general frameworks 0, [381 GEJ E7] . For these 
theories, these current radiative tests on binary pulsars are more constraining than SEP 
tests [28]. 



Nevertheless, there are two aspects that support the importance of SEP tests 
with binary pulsars. Firstly, they can be interpreted as generic, direct tests for the 
UFF of strongly self-gravitating bodies, independent of any specific gravity theory. 
Secondly, non-linear strong-field effects, like spontaneous scalarization, could be limited 
to very massive neutron stars, which until now have only been observed in wide binary 
systems like PSR J1614-2230 (M p ~ 1.97 M , P b = 8.7 days [20J) and PSR J1903+0327 
(M p ~ 1.68 M Q ,P b = 95 days [25]). In these wide systems a dipolar contribution to the 
gravitational wave damping could, even in future, be too small to be detectable, or could 
not be separable from kinematic effects [17], while the deviation from GR could still be 
measurable in an SEP test. In particular since, as outlined in Section I4~3l the e-test is 
not "contaminated" by external effects. 

6. Future prospects: Globular cluster pulsars, triple systems and mergers 

The external gravitational acceleration by the Galaxy is rather small (|g| ~ 2 x 
10 _10 ms~ 1 at the location of the Sun). In a stronger external gravitational field SEP 
violation would be proportionally more prominent, (eq. |8]). There are numerous pulsar 
binaries known to exist in globular clusters [24], where the external acceleration is 
typically two orders of magnitude larger than in the field of the Galaxy [63l [9J 126] . 
Unfortunately one cannot determine the exact location of these systems within the 
globular cluster as one does not have a good handle on the radial distance for these 
pulsars; however this can be somewhat constrained for systems with negative period 
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derivatives. Furthermore, the latter provide direct lower limits on g along the line of 
sight. If the semi-latus rectum happened to lie along the line of sight, we would then 
automatically have all we need to derive an upper limit on A p . 

An even stronger external field would arise in a hierarchical triple star system, 
where the pulsar is in a tight orbit with a weakly self-gravitating object, and this 
inner binary falls in the gravitational field of a third companion. This would resemble 
the SEP test done in the Earth-Moon-Sun system, but with a strongly self-gravitating 
object. In the globular cluster M4 there is the millisecond pulsar PSR B1620— 26 where 
the inner companion to the pulsar is a ~ 0.3 M Q WD and the outer companion is 
a ~ Jupiter mass companion [531 48J. Unfortunately, the low-mass outer companion 
yields an external acceleration g which is only a factor of a few larger than the typical 
value for the acceleration in the Galactic plane, and therefore does not lead to any 
interesting constraints. The situation would be quite different if the outer companion 
were a ~ 1 M Q star, and even better if it were also a neutron star. A future detection 
of a comparable-mass triple with a pulsar is not unlikely [54], particularly since we 
have already discovered a binary system, PSR J1903+0327, which started its life as a 
hierarchical triple system and appears to have become a binary system much later in 
its evolution [25]. Other systems with similar origins might still survive as hierarchical 
triples. 

Similarly to equation (T4j), in a hierarchical triple system one would have [TTJ ; 

A p - A c ~ a ex (a p - a c ). (22) 

According to this equation, the ideal triple system combination would be a pulsar- WD 
(a c = «o "C 1) or pulsar-black hole (a c = 0) system in the field of a more distant neutron 
star, for which a ex , like a p , could in theory be much larger than a . Consequently, in 
a hierarchical triple system not only the external gravitational acceleration felt by the 
internal binary would be much larger than for a binary pulsar falling in the field of the 
Galaxy, also the effective scalar coupling of the source of the external field could be 
significantly larger, provided the outer companion is a neutron star. 

In such a triple system the challenge lies in obtaining a sufficient number of higher 
order derivatives in the pulsar frequency, in order to constrain the orbit and mass of the 
distant companion such that any effects from the SEP violation inflicted on the inner 
orbit could be separated from "classical" orbital perturbations [SJ [33] . This depends 
on the details of the system, particularly its orbital period and the timing observations 
(precision, time span). 

Finally, if the SEP violating interaction is only of limited range, like in the massive 
Brans-Dicke theory of gravity (see [T]), then Galactic binary pulsars are insensitive 
to the corresponding violation of the UFF, since the masses that cause the external 
gravitational field are at large distances. In such a case, direct tests of strong field 
aspects of the UFF would require a hierarchical triple system. 

In about three years from now it is expected that, after their upgrade, the ground 
based gravitational wave detectors will make their first detections of gravitational waves. 
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One of the most promising sources for these detectors are inspiralling compact binaries, 
consisting of neutron stars or black holes j57j. This will not only mark the beginning 
of the era of gravitational-wave astronomy, it will also provide new possibilities to test 
gravity, including (indirect) tests of UFF for self-gravitating bodies via its theoretical 
connection to the radiative properties of a gravity theory. Although, for certain scalar- 
tensor theories (including JFBD) it has been shown that binary pulsar experiments are 
already more constraining than it is expected for the advanced LIGO/VIRGO detectors 
[621 IT5| 128] , there are many theoretical aspects where ground-based gravitational wave 
antennae will nicely complement binary pulsar experiments in the near future, notably 
theories where the gravitational interaction is partly mediated by very short range 
(~ 1 lt-s, and less) fields. 

Acknowledgements 

P.F. gratefully acknowledges the financial support by the European Research Council 
for the ERC Starting Grant BEACON under contract no. 279702. We also thank Gilles 
Esposito-Farese for useful comments and stimulating discussions and Nicolas Caballero 
for a careful review of this paper. We are grateful to Stanley Deser for bringing (29] to 
our attention. 

References 

[1] J. Alsing, E. Berti, C. M. Will and H. Zaglauer. Gravitational radiation from compact binary 
systems in the massive Brans-Dicke theory of gravity. Phys. Rev. D, 85:064041, 2012. 

[2] Z. Arzoumanian, K. Joshi, F. A. Rasio and S. E. Thorsett. Orbital Parameters of the PSR B1620- 
26 Triple System. "IAU Colloq. 160: Pulsars: Problems and Progress" , ASP Conference Series, 
Vol. 105. Edited by S. Johnston, M. A. Walker and M. Bailes. San Francisco: Astronomical 
Society of the Pacific, 1996, p. 525. 

[3] H. Bartko, et al. GRAVITY: Astrometry on the galactic center and beyond New Astronomy 
Reviews, 53:301, 2009 

[4] B. Bcrtotti, L. less and P. Tortora. A test of general relativity using radio links with the Cassini 

spacecraft. Nature, 425:374, 2003. 
[5] N. D. R. Bhat, M. Bailes and J. P. W. Verbiest. Gravitational-radiation losses from the pulsar 

white-dwarf binary PSR J1141-6545. Phys. Rev. D, 77:124017, 2008. 
[6] V. Boriakoff, R. Buccheri and F. Fauci. Discovery of a 6. 1-ms binary pulsar PSR 1953+29. Nature, 

304:417, 1983. 

[7] C. Brans and R. H. Dicke. Mach's Principle and a Relativistic Theory of Gravitation Phys. Rev., 
124:925, 1961. 

[8] D. J. Champion, S. M. Ransom, P. Lazarus, F. Camilo, C. Bassa, V. M. Kaspi, D. J. Nice, P. 

C. C. Freire, I. H. Stairs, J. van Leeuwen, B. W. Stappers, J. M. Cordes, J. W. T. Hcsscls, 

D. R. Lorimer, Z. Arzoumanian, D. C. Backer, N. D. Ramesh Bhat, S. Chatterjee, I. Cognard, 
J. S. Deneva, C.-A. Faucher-Guigure, B. M. Gaensler, J. Han, F. A. Jenet, L. Kasian, V. I. 
Kondratiev, M. Kramer, J. Lazio, M. A. McLaughlin, A. Vcnkataraman and W. Vlemmings. 
An Eccentric Binary Millisecond Pulsar in the Galactic Plane. Science, 320:1309, 2008. 

[9] N. D'Amico, A. Possenti, L. Fici, R. N. Manchester, A. G. Lyne, F. Camilo and J. Sarkissian. 
Timing of Millisecond Pulsars in NGC 6752: Evidence for a High Mass-to-Light Ratio in the 
Cluster Core. Astrophys. J. Lett., 570:89, 2002. 



Universality of free fall for strongly self-gravitating bodies 



17 



[10] T. Damour. Binary Systems as Test-beds of Gravity Theories. larXiv:0 704.0749, 2007. 
[11] T. Damour and G. Esposito-Farese. Tensor- multi-scalar theories of gravitation. Class. Quant. 
Grav., 9:2093, 1992. 

[12] T. Damour and G. Esposito-Farese. Nonperturbative strong-field effects in tensor-scalar theories 
of gravitation. Phys. Rev. Lett, 70:2220, 1993. 

[13] T. Damour and G. Esposito-Farese. Testing gravity to second post-Newtonian order: A field- 
theory approach. Phys. Rev. D, 53:5541, 1996. 

[14] T. Damour and G. Esposito-Farese. Tensor-scalar gravity and binary-pulsar experiments. 
Phys. Rev. D, 54:1474, 1996. 

[15] T. Damour and G. Esposito-Farese. Gravitational-wave versus binary-pulsar tests of strong-field 
gravity. Phys. Rev. D, 58(042001):1, 1998. 

[16] T. Damour and G. Schafer. New tests of the strong equivalence principle using binary pulsar data. 
Phys. Rev. Lett., 66:2549, 1991. 

[17] T. Damour and J. H. Taylor. On the orbital period change of the binary pulsar PSR 1913+16. 
Astrophys. J. 366:501, 1991. 

[18] T. Damour and J. H. Taylor. Strong-field tests of relativistic gravity and binary pulsars. 
Phys. Rev. D, 45, 1840, 1992. 

[19] P. B. Demorest, et al. Limits on the Stochastic Gravitational Wave Background from the North 
American Nanohertz Observatory for Gravitational Waves. \arXiv: 1201. 66^1 vl, 2012. 

[20] P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. Roberts and J. W. T. Hessels. A two-solar- 
mass neutron star measured using Shapiro delay. Nature, 467:1081, 2010. 

[21] S. Detweiler. Pulsar timing measurements and the search for gravitational waves. Astrophys. J., 
234:1100, 1979. 

[22] R. D. Ferdman, R. van Haasteren, C. G. Bassa, M. Burgay, I. Cognard, A. Corongiu, N. D'Amico, 
G Desvignes, J. W. T. Hessels, G. Janssen, A. Jessner, C. Jordan, R. Karuppusamy, E. F. Keane, 
M. Kramer, K. Lazaridis, Y. Levin, A. G. Lyne, M. Pilia, A. Possenti, M. Purver, B. Stappers, 
S. Sanidas, R Smits and G. Thcureau. The European Pulsar Timing Array: current efforts and 
a LEAP toward the future. Class. Quant. Grav., 27:084014, 2010. 

[23] M. Fierz. Uber die physikalischc Deutung der erweiterten Garvitationstheorie P. Jordans. Helv. 
Phys. Acta, 29:128, 1956. 

[24] P. C. C. Freire. Pulsars in Globular Clusters. |http:// www.naic.edu/~pfreire/GCpsr.html 

[25] P. C. C. Freire, C. G. Bassa, N. Wex, I. H. Stairs, D. J. Champion, S. M. Ransom, P. Lazarus, 
V. M. Kaspi, J. W. T. Hessels, M. Kramer, J. M. Cordcs, J. P. W. Verbiest, P. Podsiadlowski, 
D. J. Nice, J. S. Deneva, D. R. Lorimer, B. W. Stappers, M. A. McLaughlin and F. Camilo. On 
the nature and evolution of the unique binary pulsar J1903+0327. Month. Not. R. Astr. Soc., 
412:2763, 2011. 

[26] P. C. C. Freire, F. Camilo, M. Kramer, D. R. Lorimer, A. G. Lyne, R. N. Manchester, and 

N. D'Amico. Further results from the timing of the millisecond pulsars in 47 Tucanae. 

Month. Not. R. Astr. Soc., 340:1359, 2003. 
[27] P. C. C. Freire and N. Wex. The orthometric parameterisation of the Shapiro delay and an 

improved test of general relativity with binary pulsars. Month. Not. R. Astr. Soc, 409:199, 

2010. 

[28] P. C. C. Freire, N. Wex, G. Esposito-Farese, J. P. W. Verbiest, M. Bailes, B. A. Jacoby, 
M. Kramer, I. H. Stairs, J. Antoniadis and G. Janssen. The relativistic pulsar-white dwarf 
binary PSR J1738+0333. II. The most stringent test of scalar-tensor gravity. In press 
Month. Not. R. Astr. Soc, 2012, arXiv:01205.1450. 

[29] H. Goenner Some remarks on the genesis of scalar-tensor theories. arXiv:01204.3455. 

[30] M. E. Gonzalez, I. H. Stairs, R. D. Ferdman, P. C. C. Freire, D. J. Nice, P. B. Demorest, 
S. M. Ransom, M. Kramer, F. Camilo, G. Hobbs, R. N. Manchester and A. G. Lyne. 
High-precision Timing of Five Millisecond Pulsars: Space Velocities, Binary Evolution and 
Equivalence Principles. Astrophys. J., 743:102, 2011. 



Universality of free fall for strongly self- gravitating bodies 



18 



[31] J. H. Jeans. The effect of varying mass on a binary system. Month. Not. R. Astr. Soc, 85:912, 
1925. 

[32] P. Jordan. Zum gegenwartigen Stand der Diracschen kosmologischen Hypothcscn Zeitschrift fur 
Physik, 157:112, 1959. 

[33] K. J. Joshi and F. A. Rasio. Distant Companions and Planets around Millisecond Pulsars. 
Astrophys. J., 479:948, 1997. Erratum: Astrophys. J., 488:901, 1997. 

[34] V. M. Kaspi, J. H. Taylor and M. F. Ryba. High-precision timing of millisecond pulsars. 3: 
Long-term monitoring of PSRs B1855+09 and B1937+21. Astrophys. J., 428:713, 1994. 

[35] S. M. Kopeikin. Proper Motion of Binary Pulsars as a Source of Secular Variations of Orbital 
Parameters. Astrophys. J. Lett., 467:93, 1996. 

[36] S. Kopeikin and V. A. Potapov. Binary Pulsar Tests of General Relativity in the Presence of Low- 
Frequency Noise. "IAU Colloq. 177: Pulsar Astronomy — 2000 and Beyond", ASP Conference 
Series, Vol. 202. Edited by M. Kramer, N. Wex and R. Wielebinski. San Francisco: Astronomical 
Society of the Pacific, 2000, p. 117. 

[37] M. Kramer, I. H. Stairs, R. N. Manchester, M. A. McLaughlin, A. G. Lyne, R. D. Ferdman, 
M. Burgay, D. R. Lorimer, A. Posscnti, N. D'Amico, J. M. Sarkissian, G. B. Hobbs, 
J. E. Reynolds, P. C. C. Freire and F. Camilo. Tests of General Relativity from Timing the 
Double Pulsar. Science, 314:97, 2006. 

[38] K. Lazaridis, N. Wex, A. Jessner, M. Kramer, B. W. Stappcrs, G. H. Janssen, G. Desvignes, 
M. B. Purver, I. Cognard, G. Theureau, A. G. Lyne, C. A. Jordan, and J. A. Zensus. Generic 
tests of the existence of the gravitational dipole radiation and the variation of the gravitational 
constant. Month. Not. R. Astr. Soc., 400:805, 2009. 

[39] K. Nordtvedt. Equivalence Principle for Massive Bodies. I. Phenomenology. Phys. Rev., 169:1014, 
1968. 

[40] K. Nordtvedt. Equivalence Principle for Massive Bodies. II. Theory. Phys. Rev., 169:1017, 1968. 
[41] K. Nordtvedt. Testing Relativity with Laser Ranging to the Moon. Phys. Rev., 170:1186, 1968. 
[42] P.C.Peters. Gravitational Radiation and the Motion of Two Point Masses. Phys. Rev., 136:B1224, 
1964. 

[43] D. Psaltis, D. Perrodin, K. R. Dienes and I. Mocioiu. Kerr Black Holes Are Not Unique to General 

Relativity. Phys. Rev. Lett., 100:091101, 2008. 
[44] F. A. Rasio. Is there a planet in the PSR 1620-26 triple system? Astrophys. J., 427T107, 1994. 
[45] M. F. Ryba and J. H. Taylor. High-precision timing of millisecond pulsars. I - Astrometry and 

masses of the PSR 1855+09 system. Astrophys. J. 371:739, 1991. 
[46] M. V. Sazhin. Opportunities for detecting ultralong gravitational waves. Soviet Astronomy, 22:36, 

1978. 

[47] R. T. Schilizzi, P. E. F. Dewdney and T. J. W. Lazio. The square kilometre array. Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 7733:35, 2010 

[48] S. Sigurdsson and S. E. Thorsett. Update on Pulsar B1620-26 in M4: Observations, Models 
and Implications. Binary Radio Pulsars, ASP Conference Series, Vol. 328, Proceedings of 
the conference held 11-17 January, 2004, Aspen, Colorado, USA. Edited by F. A. Rasio and 
I. H. Stairs. San Francisco: Astronomical Society of the Pacific, 2005, p. 213. 

[49] I. I. Shapiro. Fourth Test of General Relativity. Phys. Rev. Lett., 13:789, 1964. 

[50] I. S. Shklovskii. Possible Causes of the Secular Increase in Pulsar Periods. Soviet Astronomy, 
13:562, 1970. 

[51] E. M. Splaver, D. J. Nice, I. H. Stairs, A. N. Lommen and D. C. Backer. Masses, Parallax and 
Relativistic Timing of the PSR J1713+0747 Binary System. Astrophys. J., 620:405, 2005. 

[52] I. H. Stairs, A. J. Faulkner, A. G. Lyne, M. Kramer, D. R. Lorimer, M. A. McLaughlin, 
R. N. Manchester, G. B. Hobbs, F. Camilo, A. Possenti, M. Burgay, N. D'Amico, P. C. C. Freire 
and P. C. Gregory. Discovery of Three Wide-Orbit Binary Pulsars: Implications for Binary 
Evolution and Equivalence Principles. Astrophys. J., 632:1060, 2005. 

[53] S. E. Thorsett, Z. Arzoumanian, F. Camilo and A. G. Lyne. The Triple Pulsar System PSR 



Universality of free fall for strongly self-gravitating bodies 



19 



B1620-26 in M4. Astrophys. J., 523:763, 1999. 
[54] M. Trenti, S. Ransom, P. Hut and D. C. Heggie. Predictions for triple stars with and without a 

pulsar in star clusters. Month. Not. R. Astr. Soc, 387:815, 2008. 
[55] W. van Straten, M. Bailes, M. Britton, S. R. Kulkarni, S. B. Anderson, R. N. Manchester and 

J. Sarkissian. A test of general relativity from the three-dimensional orbital geometry of a binary 

pulsar. Nature, 412:158, 2001. 
[56] J. P. W. Verbiest, M. Bailes, N. D. R. Bhat, S. Burke-Spolaor, D. J. Champion, W. Coles, 

G. B. Hobbs, A. W Hotan, F Jenet, J Khoo, K. J. Lee, A Lommen, R. N. Manchester, 

S. Oslowski, J. Reynolds, J. Sarkissian, W. van Straten, D. R. B. Yardley and X. P. You. 

Status update of the Parkes pulsar timing array. Class. Quant. Grav., 27:084015, 2010. 
[57] S. J. Waldman. et al. The Advanced LIGO Gravitational Wave Detector. la~rXiv:1103.2728| 2011. 
[58] J. M. Weisberg, D. J. Nice and J. H. Taylor. Timing Measurements of the Relativistic Binary 

Pulsar PSR B1913+16. Astrophys. J., 722:1030, 2010. 
[59] N. Wex. New limits on the violation of the strong equivalence principle in strong field regimes. 

Astronomy and Astrophysics, 317:976. 1997. 
[60] N. Wex. Small-eccentricity binary pulsars and relativistic gravity. "IAU Colloq. 177: Pulsar 

Astronomy — 2000 and Beyond", ASP Conference Series, Vol. 202. Edited by M. Kramer, 

N. Wex and R. Wiclcbinski. San Francisco: Astronomical Society of the Pacific, 2000, p. 113. 
[61] C. M. Will. Theory and Experiment in Gravitational Physics. Cambridge University Press, 

Cambridge, 1993. 

[62] C. M. Will. Testing scalar-tensor gravity with gravitational-wave observations of inspiralling 

compact binaries. Phys. Rev. D 50:6058, 1994. 
[63] A. Wolszczan, S. R. Kulkarni, J. Middleditch, D. C. Backer, A. S. Fruchter and R. J. Dewey. A 

110-ms pulsar, with negative period derivative, in the globular cluster M15. Nature, 337:531, 

1989. 



